Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 8 March 2013 (MN KTeX style file v2.2) 



Halo statistics in non-Gaussian cosmologies: the collapsed fraction, 
conditional mass function, and halo bias from the path-integral 
excursion set method. 
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ABSTRACT 

Characterizing the level of primordial non-Gaussianity (PNG) in the initial conditions for 
structure formation is one of the most promising ways to test inflation and differentiate 
among different scenarios. The scale-dependent imprint of PNG on the large-scale cluster- 
ing of galaxies and quasars has already been used to place significant constraints on the level 
of PNG in our observed Universe. Such measurements depend upon an accurate and robust 
theory of how PNG affects the bias of galactic halos relative to the underlying matter den- 
sity field. We improve upon previous work by employing a more general analytical method - 
the path-integral extension of the excursion set formalism - which is able to account for the 
non-Markovianity caused by PNG in the random-walk model used to identify halos in the 
initial density field. This non-Markovianity encodes information about environmental effects 
on halo formation which have so far not been taken into account in analytical bias calcula- 
tions. We compute both scale-dependent and -independent corrections to the halo bias, along 
the way presenting an expression for the conditional collapsed fraction for the first time, and 
a new expression for the conditional halo mass function. To leading order in our perturba- 
tive calculation, we recover the halo bias results of Desjacques et. al. (2011), including the 
new scale-dependent correction reported there. However, we show that the non-Markovian 
dynamics from PNG can lead to marked differences in halo bias when next-to-leading order 
terms are included. We quantify these differences here. We find that the next-to-leading or- 
der corrections suppress the amplitudes of both the scale-dependent and -independent bias 
by - 5 - 10% for massive halos with M ~ W 15 M Q /h, and ~ 30 - 40% for halos with 
M <~ 1O 14 M0 / h. The corrections appear to be more significant as the halo mass is lowered, 
though we caution that the apparently large effects we observe in the low-mass regime likely 
signal a breakdown of the perturbative approach taken here. 

Key words: cosmology: theory, large-scale structure of the Universe, inflation - galaxies: 
statistics 



1 INTRODUCTION 

The inflationary paradigm provides a robust framework for explaining key aspects of our observable Universe such as its geometric flatness, 
features of the cosmic microwave background, and the initial conditions for structure formation. Despite these remarkable successes, we still 
know very little about the physics behind it, and currently cannot distinguish among a wide variety of viable inflationa ry models. One of the 
most p romising ways to differentiate among these models is to probe the statistics of the initial density fluctuations (see lDesiacques & Seliakl 
l2010al . and references therein). While the simplest scenarios - the canonical single-field slow-roll models - predict an almost perfectly 
Gaussian distribution of initial fluctuations, more general inflationary models predict significant deviations from Gaussianity that observations 
might yet be sensitive enough to detect. There is therefore great interest in developing ways to measure primordial non-Gaussianity (PNG) 
since its detection would have profound implications for inflationary theory. 

There are presently two methods that have so far been applied with some success to place significant constraints on the level of PNG 
in our observed Universe. The first is the statistical imprint of PNG in the temperature anisotro pies of the cosm ic microwave background 
(CMB), which directly probe the initial fluctuations while they are still in the linear regime (e.g. lKomatsull2010l) . The second is the imprint 
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on the large-scale structure that develops as the initial fluctuations grow to the highly non-linear point of forming galaxies, which serve 
as tracers of the underlying matter. In general, non-linear growth can complicate the interpretation of the observed structure, both because 
density fluctuations develop non-Gaussianity even when the initial conditions are purely Gaussian, and because the theoretical predictions 
in this regime ultimately depend upon N-body simulations. On large enough scales, however, the density fluctuations filtered on these scales 
are still linear, and the problem reduces to the theory of how well galaxies trace the large-scale mass distribution - the so-called "bias." 

The prospect of probing PNG with large-scale structure improved dramatically when it was discovered that the mode coup ling effects 
of PNG induce a scale-dependent signature in the powe r spectrum of biased tracers (such as galactic halos) on large scales dDalal et al.l 
120081 ; IMatarrese & Ve rde 2008; Afshordi & Tollev 200^). The first and best studied example of this signature involves the local-quadratic 
model of PNG, in which the primordial Bardeen potential fluctuation in the matter-dominated epoch, $ng(s), is obtained from a quadratic 
transformation of the local Gaussian fluctuation field, 4>g(x), according to 

®ng(x) = 4>g(x) + /nl [<fc{x) 2 - <0|(aO>] , (1) 

where /nl is the so-called non-linearity parameter dSalopek&Bondlll990l: iKomatsu & Spergelll200lh . In this case, while the Gaussian field 
is entirely characterized by its power spectrum, P$(fc), information about higher-order correlations, for example through the bispectrum, is 
required to characterize the statistics of the non-Gaussian potential. To first order in /nl, equation (0 gives a bispectrum with the forrrQ, 

Bi° Cal (fcl,fc 2 ,fc 3 ) = 2/NL[P$(fc 1 )P*(fc2) +P*(fc 1 )P*(fc 3 ) +P*(fc 2 )P*(fc 3 )]. (2) 

Observational constraints have been placed on this form of the bispectrum by finding the range of allowed amplitudes, expressed in terms of 
/nl - For example, the CMB anisotropy measurements by the Wilkinson Microwave Anisotropy Probe seven-year data analysis (WMAP7) 
find a 95% limit of —10 < /nl < 74 dKomatsu et alj|20l"il ). The bispectrum in e quation ® is important in the phenomenology of PNG 
because a detect i on of non-zero /nl would rule out standard single-field inflation dCreminelli & Zaldarriagair2004t ISeerv & Lidsevll2005l : 
IChen et al.l2007llCheung et al.l2008h . 

In order to constrain equation (|2j from observations of large-scale structure, it is necessary to consider the expected halo bias for this 
model, i.e. the ratio of the fractional halo number density to the fractional matter density. In Fourier space, this ratio has been found to depart 
from the Gaussian expectation by a correction term, Ab(k), containing two parts: one that depends on the wavenumber (scale-dependent) 
and one that does not (scale-independen t). In the limit of s mall wavenumber, the correction approaches the fo rm Ab(k) oc /nl(&g — l)/k 2 , 
where be is the exp ected Gaussian bi as dDalal et alj2008l : lMatarrese & Verdell2008l : I Afshordi & Tolley||2008l) . Based upon this assumed k 
scale-dependence. ISlosar et al.l d2008h have already constrained /nl to be in the range —31 < /nl < 70 (95% limit) using the clustering of 
massive galaxies and quasars in the Sloan Digital Sky Survey; a result that is competitive with the WMAP7 constraints. There is naturally a 
great interest in this method as future large-scale structure surveys of e ver-increasing volume uti l izing both the galaxy power spectrum and 
bispectrum may surpass the CMB measurements in constraining /nl dScoccimarro et al1 i2004; Jeo ng & Komatsull2009l ; lNishimichi et al.l 
l2010l ; lBaldauf et alj201 lllGiannantonio et alj2012h . 

Numerical N-body methods have confirmed that the halo bias scales as k~ 2 in the small-fc limit in the local-quadratic m odel, with a 
redsh i ft dependence inversely proportional to the linear perturbation g r owth factor, D(z), as predicted by the analytica l theory dDalal et alj 
I2OO8I; IDesiacques et alj|2009l ; IPillepich et alj|20ld : iGrossi et al1l2009l ; iNishimichi et al.|[201fi ISmith & LoVerddl20lll) . However, the am- 
plitud e of the bias is still somew hat uncertain, as the N-body simulations so far produce a range of values differing at the 10 — 20% 
level dDesiacques & Seliakll2010bT) . On the othe r hand, the analytical predictions disagree significantly with results from N-body simula- 
tions in models beyond t he local-quadratic case dDesiacques & Seliakll2010cllShandera et al .11201 lUWagner & Verddl2012l) . This motivated 
IDesiacques et al.1 d201 lbl) to re-examine three analyt ical derivations of the bia s (see also lsmith et alj|2012h : 1) An approach based on the 
statistics of thresholded regions in the density field (Mat arrese & Verd"ell2008l) . 2) A peak-backgroun d split (PBS) approach based on the 
separation of uncorrelated long- and short- wavelength contribu tions to the Gaussian p erturbations dDalal et aT]|200 i lSlosaretaLll2008l; 
ISchmidt & Kamionkowskill20ld) . 3 ) A second PBS approac h dDesiacques et alj|2010l) using the conditional halo mass function, which 
they derive by an ex tension of the IPress & SchechteJ d 19741) method to non-Gaussian initial conditions (also see IMatarrese et al. I l2000l: 
ILo Verde et al.ll2008h . This approach is conceptually different from the previous one because it does not involve a separation of scales in 
the Gaussian perturbations, but instead considers through the conditional mass function how the local halo-abundance depends on the large- 
scale non-Gaussian dens ity contrast. In addition to showing that the thresholding method cannot be reconciled with N-body simulations, 
IDesiacques et al.l d2011bl) use the last two approaches to derive a ne w scale-dependent contribution to the bias that was previously overlooked 
in the literature. In a companion paper dDesiacques et al.l |20 Hal) , they showed that the new term is critical for improving the analytical 
predictions that were previously discrepant with N-body simulations. 

It is important to confirm the abo ve findings with more general methods . Here, we consider an independent analytical approach to the 
halo bias - the excursion set formalism dBond et al.ll99ll ; lLacev & Coleil993l) . Until recently, the excursion set method had been analytically 
tractable only in the case with Gaussian initial conditions and, even then, only when a sharp k-space filter w as used. In this case, the dyn am- 
ics in the excursion-set random-walk model for identifying halos in the initial density field are Markovian. iMaggiore & Riottd d2010jbHcl) 
recently showed how to extend the excursion set model to include non-Markovian dynamics by formulating it with a path integral. This break- 
through opens the door to non-Gaussian initial conditions and/or more general filter functions (for a different approach to non-Markovian 



1 The bispect rum in equation l|2) is more general than equation {T}, as it can be generated in number of different models that do not involve the latter (see 
footnote 34 in Komatsu et all201 ll . for example). It is nonetheless customary to refer to this form for the bispectrum as the local template. 
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dynamics, see lParaniape & Shethl2012l ; |Paraniape et al.l2o72llMusso & Paraniapel2012hlMusso & Shethl2012hlMusso et alj2012h. The path 



integral approach has been successfully applied in a number of contexts involving PNG {Maggiore & Riottdl2010d : lD'Amico et alj|201 id : 



Ide Simone et alJl201ld : lD'Amico et alJl2011bl : lD'Aloisio & Nat araian 20lTI: lde"simone et al.ll2011ah . However, it has not yet been used to 
derive the scale-dependent correction to the halo bias. 

In what follows, we use the path-integral excursion set method to derive expressions for the conditional collapsed fraction and con- 
ditional mass function for non-Gaussian models with general bispectra. We th en use these expressions to compute both scale-dependent 
and -independent corrections to the halo bias. In comparison with the results of iDesiacques et all ( l2011bl) . our results will come closest to 
the conclusions of their third approach described above. In the process of deriving our result, we will be able to investigate more specifi- 
cally under what circumstances their result is valid. In particular, the excursion-set approach shows through the non-Markovian dynamics 
that environmental effects on halo formation can lead to marked differences in halo bias. We will also quantify these differences here by a 
perturbative calculation, applied for illustrative purposes to the familiar local-quadratic model, as well as a second example - the so-called 
orthogonal template. 

The remainder of this paper is organized as follows. In §[2] we review statistics of the linear density field and define the bispec trum 
templates used in this work for plotting purposes. In §[3] we outline the path-integral excursion set method of Maggiore & Riottd l l2010d] bllcl). 
In § [4] we use the formalism to calculate the conditional collapsed fraction to leading-order, and also compute the next-to-leading order 
environmental corrections. We then use the collapsed fraction to obtain expressions for the conditional mass function in § [5] From the 
conditional mass function, we obtain scale-dependent and -independent linear bias parameters in § [6] Finally, we summarize our results and 
offer concluding remarks in §[7] 

When plotting our results, we use a fiducial ACDM cosmology with parameters Q m = 0.27, Qa = 0.73, Qb = 0.04 6, Hq = 
WOh km s" 1 Mpc" 1 (with h = 0.7), n s = 0.97 and a$ = 0.8 2, consistent with WMAP7 constraints dKomatsu et~aT]|201lh . We also 
employ the linear matter power spectrum of bisenstein & Huld 19991) . 



2 STATISTICS OF DENSITY FLUCTUATIONS IN MODELS WITH PRIMORDIAL NON-GAUSSIANITY 
2.1 The density contrast and its two- and three-point functions 

The excursion set model is formulated in the Lagrangian picture; it is a method for computing halo statistics from the statistics of the linearly 
extrapolated initial density fluctuations. We quantify fluctuations in the linearly extrapolated density field, p m (x, 2), with the density contrast 
smoothed on scale R about a point x, 

5 r {-k,z) = Jd 3 x' W (\x-x.'\,R)8(x.',z), (3) 

where W is a spherically symmetric filter function with characteristic scale R and the un-smoothed density contrast is S(x, z) = 
p m (x, z)/p m (z) — 1. The Fourier transform of the smoothed density contrast is related to the primordial Bardeen potential in the matter- 
dominated epoch through the cosmological Poisson equation^ 

5 R (k,z) = M R (k,z)$(k). (4) 

Here, we define 

where T(fc) is the matter transfer function normalized to unity on large scales, g(z) is the linear growth factor of the potential normalized to 
unity during the epoch of matter domination (g(0) ~ 0.76 for our fiducial cosmology), Ho is Hubble's constant, and Q m is the present-day 
matter density in units of the critical density. 

Gaussian initial density fluctuations are uniquely characterized by their variance, 

S R = a% = {5 2 R ) c = J -^M 2 R {k,z)P^{k). (6) 

More generally, higher-order correlation functions are required to characterize non-Gaussian initial fluctuations. In what follows, we consider 
non-Gaussianity that is characterized solely through the three-point function, 

{8 Rl 8 R2 5 R3 ) c = J ^^^^^^M Rl (ki, z)M R2 (k 2 , z)M Rs (k 3 , z) 

x B 9 ,(k 1 ,k2,k 3 )(2ir f5 D (k 1 + k 2 + k 3 ) cxp [-i(ki + k 2 + k 3 ) • x] , (7) 



Strictly speaking, this relation holds only in the synchronous comoving gauge jjeong et alj2012l) . 
© 0000 RAS, MNRAS 000. 000-000 
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where 8d is the Dirac delta function. We neglect the effects of all higher-order spectra. 

The excursion set method is unique among other analytical techniques for treating non-Gaussianities because it depends on 
{Sr 1 8r 2 8r 3 ) c for the full range of smoothing scales down to the Lagrangian radius of the halo. As we will see, the mixed three-point 
correlator (8r 1 8r 1 8r 2 ) c will be particularly important in our calculations; it is the source of the scale-dependence of the bias. We may 
rewrite this correlator in a way that is convenient for manifesting the scale-dependence by collapsing the delta function, rearranging, and 
relabeling the integration variables to obtain 



/d 3 k 
-j—^ Mr 2 (k, z)P^{k)AS Rl 



45 Hl P (fc)7 (2tt)3 



d 3 k 



\Mr x (k 1 ,z)MR 1 (q, z)Bg>(k, ki,q) 



(8) 



where q = — (fei + kz). The quantity in the brackets is a form factor denoted by J-"^ (fc) in the notation of lDesiacques et all J201 lbl) . We 
adopt this notation and express the mixed correlator as 



(6 Ri 5hM (x) = J -^M R2 (k,z)P^k)4S Rl F R 3) l (k). 



(2 

The advantages of rewriting {Sr 1 Sr 1 8r 2 ) c in this form will become apparent in the following sections. 



(9) 



2.2 Phenomenological templates of the primordial bispectrum 

Primordial bispectra generated by inflationary models vary considerably and can be quite complicated. For the purpose of computing the 
effects of a primordial bispectrum on large-scale structure, it is convenient to employ commonly used phenomenological templates. We will 
use the local template extensively when plotting our results, where the primordial bispectrum is given by equation (|2}. At times, we will 

(3) 

find it desirable to bring out the effects of the form factor, , in our final results. The local template is not ideal for this task because, in 
this case, J- R 3 ^ asymptotes to /nl in the low-fe limit. As a result, we will also use the orthogonal template dSenatore et al.|[2oToh . where the 
bispectrum has the form 



S ortho / j j 7 \ n i-ortho 
5 , (fcl,fe2,fc 3 ) = 6/ NL 



-3(P 4 ,(fc 1 )P $ (fc 2 ) + cyc.)-8(P 4 ,(fci)P< i ,(fc 2 )P 4 ,(fc3)) 2/3 
-3 (py 3 (fc 1 )P^ /3 (fc 2 )P*(fc 3 ) + eye 



+ ■■: 



(10) 



and /nl* 10 is a constant parameter that determines the amplitude of the bispectrum. In this case, the form factor scales as J-^ 3 ' ~ k in the 
low-fc regime. 



3 THE PATH-INTEGRAL GENERALIZATION OF THE EXCURSION SET FORMALISM TO INCLUDE NON-GAUSSIAN 
INITIAL CONDITIONS 

Here we summarize the non-Markovian extension of the excursion set formalism bv lMaggiore & Riottol (l2010allcl) . For more details on the 



original formu lation of the excursion set model, we refer the reader to the pioneering paper by iBond et al.l (119911) (also see the review of 
IZentned d2007t) and references therein). 

It is convenient to linearly extrapolate the initial density fluctuations to the present day. With this choice, the density field stays fixed 
in time, and therefore so do the variance and three-point function, while the linear over-density threshold for collapse, 5 c (z), acquires an 
additional redshift dependence, S c (z) — > S c (z) /D(z), where D(z) is normalized to unity at the present day. In the excursion set procedure, 
a filter function with characteristic scale R is centered on a fiducial point in space. The density contrast is smoothed about that point with 
some large initial scale, Rq, to obtain the smoothed density contrast, So, with corresponding variance So = ctq- In what follows, the initial 
filtering scale always corresponds to the limit, Ro oo, so that So = and So = 0. The scale of the filter function is decreased and the 
corresponding 5 and S are again calculated. As this procedure is repeated, the stochastic variable 8 executes a random walk. When 8 first 
exceeds the collapse threshold 8 C , the fiducial point is assumed to reside within a halo with mass set by the filter scale. When numerically 
evaluating our results, we will use the coordinate-space top-hat filter, where R is assumed to be the initial comoving radius of the collapsed 
over-density, so that the mass of the halo is to a good approximation given by M — 4irp m R 3 /3. The Fourier transform of the coordinate- 
space top-hat filter is 

WR(k)=3 sinikR) ~j R ; os{kR) . (id 

(kRy 

However, for simplicity, and since we are mainly interested in the effects due to non-Gaussianity, we will neglect non-Markovian correction 
terms due to the use of this filter function. In terms of the derived expressions, this is equivalent to using the sharp fc-space filter in the 
for malism. 

iMaggiore & Riott cl< l2010d lcl) showed how to formulate the above model in terms of a path integral. Consider the "trajectory" traced 
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out in (S, <5)-space. The starting point is to discretize the "time" interval [0, S] so that S n = en, where n £ {0, 1, 2, • ■ • }. The probability 
density in the space of trajectories may be written as 

W(S ; ft, ... , S n ; S n ) = (5 D (S(S 1 ) -St)... S D (5(S n ) - S n )). (12) 
The integral representation of the Dirac delta function, 



8d(S)= I ^e, (13) 



is then used to write 



W(8 ;S 1 ,...,S n ;S n ) = / PAexp Aift + ^f- . . . A^ . . . Xi p {8( S ii) • • ■ - (14) 

V i=l p=2 ^' i 1= l i p =l / 

where we have used the notation, 

From here on, we will use the notation & = 8(Si). The probability density for a trajectory to obtain a value <5„ at the independent coordinate 
S n , without having crossed 5 C , is obtained by integrating W over the space of possible trajectories, 

U e (5 ;6 n ;S n )= f " dS ± . . . d<5 n _! W(8 ; Si, . . . , S n ; S n ). (16) 

Note the dependence of II e on the connected correlators for scales spanning the "time" interval. Through this dependence, the path-integral 
excursion set formalism contains additional information about the relationship between halo formation and the environment compared to the 
classic Press-Schechter approach and its simple extensions to the case with PNG. 

The upper limits on the integrals in dl6l limit trajectories to stay below the collapse barrier. In the excursion set model, the fraction of 
mass contained within halos with masses greater than M n is equal to the fraction of trajectories that cross S c before S n (M n ). This fraction 
is calculated by integrating d 1 61 > over all possible 8 n to obtain the fraction of trajectories that reach S n without having crossed S c , and then 
taking the complement, 

/coii(Sy = 1 - [ C dS n n E (<5 ; 8 n ; S n ). (17) 

J — oo 

After taking the continuum limit, e — > 0, the mass function follows from 

dn p m d/coii dS 



dM M dS dM 



(18) 



So far the discussion has been completely general. Let us now consider the evaluation of the above expressions in specific cases of 
interest. First, we take the case of Gaussian initial conditions, where the connected correlators with N ^ 3 vanish, and we are left with 

W(5 ;...,S n ;S n )= f VX exp ( i£ Xfa - i £ AiAj<<Wi}c ) ■ (19) 

J V »=1 i,3=l / 

The situation simplifies even more if we use the shaip k-space filter function, for which {SiSj) c — min(S'i, Sj). In this case let us define 

W sm (8 ;...,5 n ;S n ) = /©A expJi^A.^-i ^ A i A,min(S i ,S j ) j. (20) 
iMaggiore & Riottol d2010al) have shown that the con'esponding probability density, 



(21) 



nf»(<Mn;Sn)= \ C dS 1 ...dS n - 1 W gn \S ;...,S n ;S n ) y 

J —oo 

yields the standard excursion set result upon taking the continuum limit (e —> 0), 



U^ (S ;S n ;S n ) = 



\/2nS n 



&l\ ( (25 C -S 



cxp I 7-5— I — exp 



25V. / V 2*3. 



(22) 
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In equations J20b . ( 121b and {22}, the superscripts "g" and "m" represent Gaussian and Markovian respectively; the latter meaning that the 
correspondin g random walks are Markov ian processes - a property exhibited only in the case of Gaussian initial conditions and the sharp 
k-space filter. iMaggiore & Riottd 11201 Oaf) have also shown how to handle non-Markovian effects from the coordinate-space top-h at filter in 
an ap proximate way by perturbing W over the Markovian expression, though we do not consider such effects here (see however llVla et al.l 

Eq£i|>. 

Let us now turn our attention to the case with PNG - the focus of this work. As stated in the last section, we neglect the effects of the 
higher-order connected N-point functions, where TV > 3 (i.e. we consider only the primordial bispectrum from PNG in our calculations) so 
that 

W(5 ;...,6 n ;S n )= f VX exp (i ^ AA - ^ J2 ^<^) c ) x exp [ J2 XiXjXk (SiSjS^A . (23) 

J V i=l *,3'=1 / \ i,j,k=l ) 

From here on we will use the notation <9; = d/dSi. Expanding the second exponential in {23} under the assumption of small {Si8jSk) c and 
using the property, A exp(iX8i) = — idi exp(iX8i), equation {23} becomes 

1 

W(S ;...,5 n ;S n )^W gln (5 ;...,8 n ;S n )~- V (<Wi<5*}c 8,8, W^*,; . . . , 5 n ; S n ). (24) 

6 f—' 

This expansion^ fMaggiore & Riotto 2010c]) forms the backbone of our perturbative calculation of the collapsed fraction described in the 
next section. 



4 THE CONDITIONAL COLLAPSED FRACTION 

Consider a large, spherical region of the Universe with Lagrangian radius R = Ri, containing mass Mi, and with a smoothed linear density 
contrast, Si. The conditional collapsed fraction, fcoil(S a \5i, Si), is the fraction of mass in halos with masses between some lower threshold, 
M s (corresponding to smoothing scale R — R s ), and Mi. It is derived within the excursion set formalism from the conditional probability 
for a trajectory pass ing through (Si, Si) to reach (S s , S„), without ever having crossed th e barrier S c . This probability density is represented 
in discrete form by {Ma et alj201 ll ; |P' Aloisio & NataraiarfeOl ll ; lde Simone et alj|201 lal) 

p,, n ir q \ P(Si,Si,S 3 ,S s ) X!L dgi . . . d6^i dS l+1 . . . dt^rW (Sp; S u . . . ,S a - S a ) 

f(0 s , Os\0l, Ol) — =-z-z — x = -j . (23) 

PtfhSl) J 6 J oo dS 1 ...dSi- 1 W(S ;S 1 ,...,Si;Si) 

On the right hand side of {25}, "s" and "I" are integers given by s — S s /e and / = Si/e respectively and I < s. The conditional collapsed 
fraction corresponds to the fraction of excursion-set trajectories that first cross 5 C between the scales Si and S s , which can be obtained by 
integrating equation {25} over d8 s and taking its complement, 

fcon(S.\6i,Si) = 1- / " d6 s P(S s ,S s \8i,Si). (26) 



We start the evaluation of equation {25} by expanding W(5o; . . . , S s ; S s ) in the numerator and keeping terms up to first order in the three- 
point correlator, 



6 



/[ s s (—i) 3 

DAexp I iJ2 x i S i ~ g x i x A S i S j)c + A * A 3' A * ( S i S i S k)c 

w W gln (6 ;...,6 s ;S s )-l Y (6i6 i 8 k ) B a i d i d h W* a (6o;...,6.;S a ), (27) 
6 f— J 

where the first terms on the right-hand side correspond to the Gaussian contribution and the non-Gaussian collection respectively. Similarly, 
in the denominator of {25}- we have 

i ; 

W(5 ;...,Si;S,)^W Bln (6 ;...,Si;Si)- - (S^S k ) c 8 l d J d k W gln (S ; . . . ,8 r , Si). (28) 

In the numerator, we perform a separation between the density contrast smoothed on the halo (R a ) and environmental (Ri) scales. This 
is achieved by breaking up the sum in the second term of {27}, 

3 iD'Amico et alj fcpl lal) point out in their mass function calculation that the second term on the right-hand side of equation {24} yields a term that goes as 
<5 3 (<5 3 ) c /cr 6 , which can be of order of unity for very large masses (see their Figure 2), indicating that the expansion breaks down in that regime. 
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Figure 1. The three-point connected correlators, (<5f } c and (<5g<5;) c , which appear in the leading-order conditional collapsed fraction for the local (left) and 
orthogonal (right) templates. We assume a fixed Mi = 10 16 M©/fr. (R L = 45.3 Mpc), corresponding to 5; = 0.07. 



s I Is 

^ {Si6j5 k )c didjdk = ^2 {^3$k)c djdjdk + 3 ^ {$i$j8k)c didjdk (29) 

i,j,k=l i,j,k=l i,3=l fc=i+l 

Is s 

+ 3^ 22 {5i5j 8k) c didjdk + 2Z {Si8j8 k ) c didjdk- 

i=l j,k=l+l i,j,k=l+l 

The remaining details of the collapsed fraction calculation are reserved for the appendix due to their technical nature. An important detail, 
however, is that the triple sums in equation ( 129) become triple integrals over Si, Sj, and Sk in the continuum limit, which ar e not ana- 
lytically tractable. However, w e c an move forward analyti cally with an extension of the approach of Maggiore & Riottol d20 1 Och (also see 
ID' Aloisio & Nataraianl J20TH) and lde Simone et all d201 lah ). We employ a Taylor expansion of (SiSjdk)c about the point (S a , Sb, S c ), 



(6 i 8j5 k ) c = V ( iy+«+r _ _ _ 

z — ' p\q\r\ 

p,q,r-0 e * 



QP Qq QT 

dsfds]dsl {5lS3Sk)c 



(30) 

Si=S a ,Sj=Sb,Sk=S c 

For illustrative purposes, in Figure[T]we show the three-point connected correlators, {8^} c and (3^8i) c , for Mi = 10 16 M© /h, corresponding 
to Si — 0.07 in our fiducial cosmology. In the remainder of this work we will consider the terms in ( I30t originating from the p = q = r = 
and p + q + r = 1 contributions. In keeping with the terminology of iMaggiore & Riottol yOlOcJ), we will refer to the p = q = r = terms 
as "leading-order," while the p + q + r = 1 terms will be called "next-to-leading order." We again emphasize that all results in this paper are 
to first order in the three-point correlator. 

4.1 The leading-order terms 

To obtain the leading-order terms in equation ( 129) , we expand each of the correlators about the endpoints of the sum in which they appear, 
and keep only the constant terms with p = q — r = 0, so that ( 129) becomes 

s I Is 

]T (8iSj5 k ) c didjd k « {6f) c didjdk + 3(6 S 6?) C ^ didjd h (31) 

i,j,k = l i,j,k — l i,3 = 1 k — l-\-l 

I s s 

+ 3(5 S 2 ^) C ^ J2 9idjd k + {5 3 s ) c J2 d ^ dk - 

i=l j,k=l+l i,j,k=l+l 

The remaining steps of the calculation are detailed in Appendix lAll The final expression for the conditional collapsed fraction up to leading 
order is 
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Figure 2. Left panel: the conditional collapsed fraction up to leading-order at z = for the local template of PNG as a function the lower mass-threshold, 
Ms, for fixed environmental mass M;. We assume an environmental density fluctuation of mass Mi = fO 16 M //i (Ri = 45.3 Mpc), with peak-height 
v t = &l/V~Sl = 1- The bottom panel shows its ratio to the Gaussian case. Right panel: same as the left panel, but as a function of the environmental 
peak-height, for fixed M 3 = 1O 13 M /h. 



t _ fg m I ft 1 ) 
J coll — J coll "r" J coll ' 



(32) 



where 



i) Afs c -s l i \ df!Z , b 



f (i) _ 

J coll — 



3 \ S s — Si 5 C — Si J dSs Si 
+ C- 



5c — (5c — Si) coth f 



51 - 5 c 5i 



8 f gm 

"^coll 

dS s 



>S?-Si-2(5 2 c-5c5i) 



Sf(5 c -5 t 

and the connected correlators evaluated at the two scales enter through 



coth 



51 - 5 c 5 t 
Si 



dS s ' 



(33) 



A = A(Si,S s ) = {5'i)c - (5?}c + 3 {5?5 3 )c - 3 (5i5 2 s ) c 



B = B(S,,Ss) = (5f)c + (5i5 2 s )c - 2 {SfS a ) c 



(34) 
(35) 



C = C(S h S a ) = (5?5 a )c-(5f) c 



(36) 



Here, the super-script in the second term on the right-hand side of equation l |32t denotes the term that is leading-order in the expansion l |30| l, 
while 



/c g oTi = erfc 



5r — 5l 



V2[S S - Si] 



(37) 



and 



"J coll 
8Ss 



{6c -Si) 



^(S s - Si)*/ 2 



exp 



(5c - 5 t ) 2 



2(S a - Si) 



(38) 



are the standard excursion set expressions for the conditional collapsed fraction and first-crossing rate in the case of Gaussian initial conditions 
and Markovian random walks. Note that l |32t takes the form of the Gaussian and Markovian result plus correction terms proportional to the 
three-point connected correlators. 
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Figure 3. The effects of the next-to-leading order (NLO) corrections to the conditional collapsed fraction. All parameters are the same as the left panel of 
Figureff] The top panel shows the ratios of both the leading-order (thin, dashed and dotted) and next-to-leading order (thick, solid and dot-dashed) cases to the 
Gaussian and Markovian case, while the bottom panel shows the fractional change from adding the next-to-leading order corrections. 



Figure[2]shows the conditional collapsed fraction at z = to leading order for the local template. The top panel on the left shows / co n 
as a function of M s for a typical large-scale mass fluctuation with Mi = 10 16 Mq /h. Here, what we mean by a "typical" fluctuation is that 
the environmental peak height, Ui = 8i/\fSl is unity. For Mi = 10 16 Mq //i, this translates to 5i = 0.27 in our fiducial cosmology. PNG 
has a larger impact on the conditional collapsed fraction when only larger halos are included in the latter (case with larger M s ). PNG has a 
smaller effect when the conditional collapsed fraction is dominated by less massive halos (case with smaller M 3 ), whose abundances are not 
as strongly affected by PNG. The plot on the right of Figure|2]shows / coU for a fixed M s = 1O 13 M /h and Mi = 1O 16 M /h as a function 
of the environmental peak height, i>i. It shows that PNG has a larger impact on the conditional collapsed fraction within higher significance 
(i.e. rarer) peaks. 



4.2 The next-to-leading order terms 



The next-to-leading order terms in the conditional collapse d fraction are obtained from the p + q + r — 1 terms in equation yOJ- We adopt 
a short-hand notation similar to Mag giore & Riottol ( 1201 Pel) for the derivatives of the three-point correlators, 



G 



a,b,c 



QP Q1 Qr 



Si — SajSj—Sjj^k — Sc 

The details of the calculation are given in Appendix I A2I The next-to-leading order terms in the conditional collapsed fraction are 



(39) 



f(2) _ O ^(i,o,o) r f(i,o,Q)\ ~ Si){28 c - Si) 

./coll — Z \'~'l,l,l '-'1,1, s I g 

2G ^^_Si Sc{5c 



(r (i,o,o) 9r ,(i,o,o)s (8c - 8i) 2 

Z ^l,l,s ) g _ g 



,0(1,0,0) (8 



fi -Si 



1 — coth 



+ 



(0,0,1) 



2G 



G 



(0,1,0) 
l.s.s 



G 



(0,1,0) 
l.s.s 



(0,0,1) <5c(5c — 81) 2 



S, 



S c — 8 c 8i 
Si 

— coth 



(S 3 - Si) df t 



coll 



(5, 



8 c Si 



Si 



dSs 



rgm 
J coll 



G 



(0,0,1) 
1,1, a 



G 



(1,0,0) 



G 



G 



V2?r(B - Si) 
-0(0,0,1) 

(0,1,0) Ol ^I.i.s 



Si 



2Si 



+ 



3Si -8 C + 



(0,0,1) 
1,1, s 



6 -Si 
2tt 



1 - 



2n(S s - Sif df!™ 
(8 C - Si) dS a 



(S c -Si)Si' 
Si 



fgm 

J coir 



(40) 
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Figure [3] shows the effect of adding the next-to-leading order terms in the local template. As in Figure [2] the environmental mass and peak 
height are set to Mi = 1O 16 M0//i and t/; = 1 respectively. The top panel shows the ratios of both + f^~L (thin, dashed and dotted) 
and + / c ^u + /^jj (thick, solid and dot-dashed) to the Gaussian and Markovian collapsed fraction, /^n- The bottom panel shows the 
fractional change from adding 2 1, indicating in this example that the next-to-leading order corrections act to suppress the leading-order 
result by ~ 2% for /nl = 100 and to enhance by the same factor for /nl = —100. 



5 THE CONDITIONAL MASS FUNCTION 

The conditional mass function is an important ingredient for deriving the halo bias in the next section. We therefore devote this section 
towards obtaining leading-order and next-to-leading order expressions for it. The conditional mass function may be written in terms of the 
first-crossing rate, 



dn c _ p m d/coii dS s 
dm m dSs dm 



(41) 



The middle factor on the right-hand side requires taking a derivative with respect to S a of the collapsed fraction that we obtained in the last 
section. Let us now consider the leading-order contribution, which may be written as 



° J roll 

dS a 



o rgm 
U J coll 

dS 3 



( OA ,\ f dB „\ . ( dC 



(42) 



where we define the following auxiliary functions to compactly express the conditional mass function, as well as greatly simplify the bias 
calculation below: 



(5c 



Si? 



2(S S -Si) 2 2(S S -Si) (S c - S t ) 2 - (S a - Si) 
(5 c -5,) 2 -3(S B -Si) 



2(S a - St) 2 

. (5c - 8j) 2 - (S B - Sj) 
2(S a - Si) 2 

1 f 5 C - 5, 1 



3 I S s — Si 5 C — 5i 
5c-(5c-5i)coth^^^ 



(43) 



(44) 



(45) 



(46) 



(47) 



S, 



Sf(5c-5i 



5i — Si — 2(5 C — 5 c 5i) 



coth 



5 c 5i 



S, 



- 1 



(48) 



At first glance, such auxiliary definitions might seem cumbersome, but the payoff is apparent for those who follow the details of the bias 
calculation given in the appendix. Using this notation, the conditional mass function up to leading-order is 



dn c 
dm 



dn c 
dm 



where 



(49) 



dn c \ 
dm / 



m dS a 



dS s 



dm 



(5 C 



■ exp 



(5c 



2(S S - Si) 



dS 3 



dm 



(50) 



m y/2n(S s - Si)*/ 2 

is the conditional Press-Schechter mass function dLacev&Cole1ll993h. Aside from o ur restriction to the spherical collapse model, we note 
that our expression differs from the moving-barrier result s of de Simone et al. U2011ah . since we include terms involving {5 S 5 2 ) C and (5f ) c , 
which are important unless Si is very small. The upper panel of Figure |4]shows equation J49t at z = for fixed Mi = 1O 16 M0 /h, vi — 1, 
and /nl = ±100, while the bottom panel shows their ratio to the Gaussian case. Before moving on to next-to-leading order contributions, 
we note that the unconditional non-Gaussian mass function can be obtaine d from equation l!49t b y taking the limit Si — > and 5"; — > 0. 
In this limit, equation l |49l > recovers the mass function originally derived bv lLo Verde et al.l d2008l) using the Edgeworth expansion and the 
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Figure 4. The conditional mass function up to leading-order at z = for the local template of PNG. The bottom panel shows the ratio of the non-Gaussian 
cases to the Gaussian case. 



Press-Schechter approach (i.e. integrating the probability density function for S), though the apparent correspondence of their result with the 
excursion set result is not triviaQ 

A full calculation of the next-to-leading order terms would require differentiating equation d40t with respect to S B . Fortunately, we do 
not need to perform this cumbersome task for the bias calculation below. There, we may focus on the large-scale environment regime, which 
is much simpler in practice because many of the above terms disappear in that limit. For this reason, we do not consider the full next- to- 
leading order corrections to the conditional mass function here, but instead consider the case where \Si\ <C 1 and Si <C 1. In this regime, we 
keep only terms in the leading-order expression involving correlators that are linear in Si, so that A -> (S 3 S ) C - 3(Si5 2 3 )c, B -> (5iS 2 a ) c , and 
C — > 0. The next-to-leading order corrections to the first-crossing rate immediately simplify because of the 1 — coth [(S 2 — 8 c 5i )/Si] factors 
appearing in l !40t , which are zero in the limit of Si — > 0, and suppress any terms they multiply. Of the remaining terms, we neglect those 
involving G; /S; since they are several orders of magnitude smaller than the Gi^i ' and Gj^J 1 , terms when \Si\ <C 1 and Si <C 1. 
With these simplifications, the conditional first-crossing rate up to next-to-leading order in the large-scale environment regime (\Si\ <C 1, 
Si < 1) is 



where 



and 



a/con _ 0/S + df^ + df& 



,(2) 



dS s 



dS s 



dS s 



dS s 



- 3 + a(S a ) c - 3a(6 s Si) c x + ' + /3 (8 3 6i) c %j) 



dS s dS B 



V dSs 



dS B 



V dSs 



(51) 



(52) 



U J coll 



2G 



((].: 



(~<\L,<J,U) \ »/v. J coll _|_ n 

ds s " v-" 1 ."." ~ ~ ~ ) 2 ds 3 \ ds a ds s 

For convenience later on, we have defined another auxiliary function, 



M) _ mi,o,o)\ 3X dfZ ^ f-ggg _ aC&T \ dfZ _ , 0>1 ,o) Si dfZ _ dG %T s. 



dS s 



G 



l rgm 



Si dSs 



dS s 



31 



coll- 



(53) 



/' : 



S s — Si 
S c — Si 



(54) 



4 See Maggiore & Riotto 1 20 1 0c) for a detailed comparison of the Press-Schechter and excursion set approaches in the context of PNG. 
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and have used the fact that 



d f.dfSA _3xdfS: (55) 



8S S V dSs J 2 8S S 

to obtain equation d53t . The next-to-leading order conditional mass function in the large-scale environment limit then follows from d41t and 
( 15 U . We will use equations d52t and d53t to derive the halo bias in the next section. 



6 THE LINEAR HALO BIAS 

We now use the results of the preceding sections to obtain the linear halo bias up to leading and next- to-leading orders in d 30b - The Lagrangian 
overdensity of halos is 



dn c f dn\ 1 df con ( d/ co ii,o \ 1 . 

Sh= d^{a^) ~ 1 = ^s7{^Sr) - 1 ' (56) 

where the "0" subscript denotes the unconditional first-crossing rate, obtained from evaluating df co \\/dS 3 at Si — and 5i = 0. Expanding 
equation d56t and keeping to linear order in the three-point function yields 



a rgm 
"./coll 



/ 8 f sm 

j "Jcoll.O 



dS 3 V 9S S 



+ 



8 f gm 

°Vcoll,( 

8S S 

QJ-gm 

8S S 



Of™ (df?°? h0 

dSs V dS s 



gm 
coll,0 



8Ss V dSs 



2 df ( c lo dfZ 
8S S 8S a 

" 2 8f^ ^)f gm 
8Ss 8Ss ' 



(57) 



The scale-independent c ontribution comes from terms containing (5j) c , while a scale-dependent contribution originates from {8^Si) c (i.e. 
iDesiacques et al.l2011bl) . We now treat each case separately. 



6.1 Scale-independent contribution 

First consider the terms up to leading-order, which are given by the first bracketed term and those denoted with superscript of (1) on the 
right-hand-side of d57t . The first bracketed term yields the old excursion set expression for the linear bias with Gaussian initial conditions. 
This may be seen by expanding it to first order in Si, and taking the limit as Si — > 0, 



lim 

Si-vo 



°Vcoll,0 

8S S 



85, 8S S 



(58) 



The second term in d57t gives both scale-dependent and -independent contributions, since it contains both {Sg8i) c and {8g) c , while the third 
term gives only scale-independent contributions. Extracting only the terms yielding scale-independent contributions, we have 



<r(i,l) 



8S S 



( dfSlo 
V dS s 



d(5p c 
8S S 



(X ~ Xo) + {$s)c (ax - aoXo) 



(59) 



where the superscript (i, 1) denotes leading-order scale-independent terms, and "0" subscripts on the a and x again indicate evaluation at 
Si — and 8i = 0. Following the same procedure used above (see Appendix IB 1 . 1 1 for details), we obtain the leading-order non-Gaussian 
scale-independent bias, 



Ab 



(i) 



,(3) 



-sr. 



8 In S s 



+ 



1 



381 
SI 



2 



(60) 



,(3) 



where is the skewness, 5i 3 ' = {8^ )c/Sl . Note that this result is identical to equation (115) of IDesiacques et all d201 lbl ) for the iV = 3 
case. We have recovered their result from the leading-order excursion set method, where the connected correlators appear in integrals over 
scales down to the halo scale, but have to a first approximation been treated as constant functions of the smoothing scale. 

Let us now consider the next-to-leading order corrections to d60t . Referring back to equations d53) and d57t , the scale-independent 



contributions come from sub-terms within the last two terms of d57| > involving either G 



(1.0,0) 



oi 8G 



(1.0,0) 



/8Ss 



?(i,2) 



J coll 



"J coll, 



3G 



(1,0,0) 



h 8S S V dS a 

Following the usual procedure (see Appendix IB 1 .21 , we obtain 



(Xo 



9G (1 '°' 0) 



(61) 
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Figure 5. Left panel: contributions to the scale-independent bias correction at z = in the local template of PNG as a function of halo mass. The solid, dashed, 
and dotted lines show the leading-order, next-to-leading order, and sum respectively. Right panel: ratio of the next-to-leading to leading order contributions to 
the scale-independent bias correction. 



o(3) 



(3) 



d In S. 
d In S s 



7_ l_ 

'6? + S s 



2Si d z S 



2 C (3) 



2 2 

5f + sT esi 3) ds * 



(62) 



Note that in addition to corrections to the existing leading-order terms, there is a new contribution involving a second derivative of the 
skewness. Finally, adding this result to d60t gives 



Ah 
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-S a 



a In 5, 



(3) 



dlnS s 



9_ 1_ 

S 2 + S s 



_3_ 3£| 

+ 81 + S* 



Ss 5?S. 



2si d 2 s 



2c(3) 



( 3 ) os 2 



(63) 



for the scale-independent bias up to next-to-leading order. 

Figure [5] shows leading and next-to-leading order contributions to the scale-independent bias at z = 0, as well as their sum for the 
local template. Clearly, the next-to-leading order corrections are sub-dominant, relative to leading-order in this case of larger halo masses, 
suppressing the amplitude by ~ 5% at M ~ 10 Mq/H. The corrections appear to become relatively more important as the mass is 
decreased. As we discuss below, this likely signals a breakdown of the expansion, 130t , 



6.2 Scale-dependent contribution 

The scale-dependent contributions to the bias come from the second and fourth terms in ( 157) . However, equati on (157) is written in coordi nate 
space, whereas the scale-dependence is manifested in Fourier space. We employ a convenient trick used by iDesiacques et alj d201 lbh for 
the conversion. Their strategy is to take the cross-correlation, (5h8i), between the halo number over-density and the large-scale smoothed 
density contrast, and rearrange the equations to pick off the bias in Fourier space, b(k). Although the excursion set model does not formally 
incorporate a p eak constraint, we emp l oy a li near halo bias relation that follows from the correlations between density peaks, 5h (k) = 
b(k)8 s (k) (see IDesiacques et alj|201ol, 1201 lbh . They showed in the case with Gaussian initial conditions that the peak constraint yields 
scale-dependent bias corrections proportional to k 2 and higher powers of k. We can therefore neglect the effects of the peak constraint since 
we are focused on the small-fc regime. Using the above bias relation, we write the left-hand side of l |57t as 

{ShS^ = J A^b(k)M s (k)Mi(k)P4,(k). (64) 

If we can now rewrite the appropriate terms on the right-hand side of l |57t in a similar way, then we can simply read off the bias coefficient 
in Fourier space. 
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Figure 6. Effects of the next-to-leading order corrections to the scale-dependent halo bias at z = in the local template of PNG. Left panel: the scale- 
dependent bias to leading order (light and thin) and next-to-leading order (dark and thick) for fixed halo masses of 10 15 M Q /h (solid) and 10 14 M Q //i 
(dashed). Right panel: fractional next-to-leading order correction as a function of halo mass for several different wavenumbers. 



Referring again to equations 1 1521 and 1 1571 . the leading-order scale-dependent contributions come from 



5 p 



dS a 



dS s 



SS S 



r,d{5 a Si} c /j2i , 
3_ Tjg — X ~ 3{dJ t ) c a X 



(65) 



The remaining steps are detailed for Appendix (TB2.ll . In summary, we expand each of the terms to first order in 8i and take the cross 
correlation with <5;. In doing so, the zeroth order terms in each expansion vanish since (5;) = 0. The cross-correlations with the linear terms 
yield factors of Si. We then take the limit as Si — > of the coefficients of the correlators and find that there is a contribution from the first 
two terms on the right-hand side of ((65}, while the last two terms give zero. Adding up the contributions from the first two terms, we find that 
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d 3 fc 

(2n? 



/^\k)Ms{k)Mi{k)P^k) = 
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2S S 
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S s 



(66) 



where have expressed the cross-correlation of 8i with the left-hand side of ([65} using equation ( 1641 ). We can now substitute (8'i8i) c from 
equation ([9} to obtain 



(2tt) 3 
which implies 



k Ab ( ^(k)M s (k)Mi(k)P4k) 
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(2tt) 3 
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S s S c 
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(67) 



A6«(fc) = 2 



( <^ _1\ 2 91nJ-j 3) (fc) 
S s 8 C J d In Ss 



■7T J (fc) 
Ms(k) 



(68) 



for equation d67t to hold for all k. This expression is identical to equation (118) of iDesiacques et all d201 lbl) for the N = 3 case. In the 
local template, the form f actor (J r j 3 ' > ) approaches /nl for small fc, and the first term within the brackets yields the well-known 1/fc 2 form for 
the scale-dependent bias jDalal et al 1 l2008l : I Matar rese & Verddl20O8t lAfshordi & Tolievll2008l) . We have also recovered the additional term 
involving dlnJ r i 3 \k)/dlnS s , which was only recently pointed out bv lDesiacques et all d201 lbl) . Although this term is s mall in the low-fc 
limit in the local template, and therefore of minimal consequence in that case, it is important for other templates as shown in lDesiacques et ail 
d2011al) . 

From equations d53t and d57t , the next-to-leading order corrections come from 
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Figure 7. Same as Figure[6] but for the orthogonal template of PNG. 
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Applying the same procedure as above, the terms in the first set of brackets vanish, while the last two terms give (see Appendix |B2.2| l 
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(70) 



We add this result to equation ((68} to obtain perhaps the main result of this paper; the scale-dependent bias up to next-to-leading order in 



dln^ (3) (fc) C, a 2 ^i 3) (fc) l .Fi 3) (fc) 

(71) 

Figure|6]shows the scale-dependent bias in the local template to leading order (light and thin) and next-to-leading order (dark and thick) 
for fixed halo masses of 1O 15 M0 /h (solid) and 10 14 Mq fh (dashed). The right panel show the fractional next-to-leading order correction as a 
function of halo mass for several different wavenumbers. In the local template, the correction is mainly due to the first term in ( |71t . The other 
terms play an insignificant role for the wavenumbers considered. In Figure|7] we reproduce Figure|6] but for the orthogonal template, where 
the other corrections in ( 171) involving derivatives of the form factor make a significant contribution. The next-to-leading order corrections 
suppress the amplitude of the scale-dependent bias by ~ 5 — 10% for halo masses of M ~ 10 1 Mg//i, and appear to get as large as 
~ 20 - 40% at M ~ 10 14 M Q //i. 

Interestingly, the first term in J7H contains a correction that can make a significant difference, even in the low-fc limit of the local 
template, particularly as the halo-mass is decreased. In this case, the above expression reduces to 



Ab d (k) = 2 15, 



+ 



2 /TO,o 
Ss [ dS a 



rgm 
Jcoll.Q 



A , /,x _ (5c 2 \ 2f NL 8 c ( E 1 \ 

Abd{k) " JMk) U - T c ) = JUk) V^'-TJ- (72) 

When the mass (and therefore &q) is decreased, the next-to-leading order correction, 1/5 C , leads to an increasingly suppressed amplitude 
compared to the case without correction, as illustrated Figures|6] That a supposedly sub-leading correction could make such a big difference 
in the low-mass regime is likely a sign that the evaluation of the path-integrals through the expansion in equation d30t may not be accurate in 
this regime fsee lMaggiore & Riotto 2010c). We therefore suggest that the scope of equation d71t should not be extended to smaller masses; 
rather, it should be limited to rare halos. Despite this possible limitation, equation l |71| > illustrates, for the first time, next-to-leading order 
effects on the scale-dependent bias from PNG within the excursion set formalism. Finally, we note that the correction in equation {72} is equal 
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to the merger bias term analytically derived by Isiosar et"al1 d2008l) . However, we stress that our correction has nothing to do with mergers, 
but instead comes from environmental effects encoded in the non-Markovianity from PNG in the excursion set random walks. 



7 SUMMARY AND DISCUSSION 

We have used the path-integral formulation of the excursion set method with constant barrier and sharp fc-space filter to calculate the condi- 
tional collapsed fraction, conditional mass function, and linear halo bias, including both scale-dependent and -independent contributions, in 
the case of non-Gaussian initial conditions with general bispectra. Our main results are summarized as follows: 

• Equations H2\ and d33 1 give the collapsed fraction to leading order in the expansion of the three-point connected correlator in equation 
d30t . In addition to being useful for other applications, this expression serves as a starting point for deriving the leading-order conditional 
mass function and linear halo bias. 

• The next-to-leading order term of the collapsed fraction is given by equation J40t . When combined with the leading-order contribution, 
this term more accurately accounts for the hierarchy of three-point correlators appearing in the excursion-set random-walk distribution 
function, equation M6\ . 

• O ur expression for the leading-order conditional mass function is equation ((49}. We note that this equation recovers the lLo Verde et alj 
d2008h mass function in the limit where the environmental density filtering radius is arbitrarily large. Next-to-leading order corrections in the 
large-scale environment limit of the conditional mass function can be obtained from equations ((52} and ((53}. 

• The leading-order scale-independent and -dependent linear bias parameters ar e given by equation (160) a nd ( 168) respectively. The path- 
integral excursion set method to leading order recovers equations (115) and (118) o flDesiacques et ail fcOllbl) for N = 3. Notably, we repro- 
duce th e additional scale-dependent term involving a derivative of the form factor, which was only recently pointed out bv lDesiacques et all 
d2011bh . While this term is negligible for the local template in the large-scale regime, it is importa nt in other templates, wh ere neglecting it 
significantly degrades the accuracy when compared to the bias measured from N-body simulations dDesiacques et aklboilah . 

• The scale-independent and -dependent bias up to next-to-leading order in ( 130) are given by equations ( 163) and ( 171) respectively. For 
cluster-halo masses (M ~ lO 15 M0/fe), the next-to-leading order terms suppress the amplitude of both contributions to the bias by ~ 
5 — 10%. In both cases, the relative effects of the next-to-leading order corrections grow as the mass is decreased. However, we caution that 
such large effects in the low-mass regime (see Figures [6] and [7] for example) may signal the breakdown of the expansion in ( 130t . We note 
that even the amplitude of the scale-dependent bias of the local template in the low-fc limit is modified in the excursion set prediction, thus 
altering a well-known result (see equation (172)). 

For the bispectrum templates tested in this work, the next-to-leading order terms become similar to the leading-order terms for masses 
as large as M ~ 5 x 1O 13 M0 /h (for which y/S/8 c = 0.62). This calls into question the validity of the expansion OOb at those masses 
and below. On the other hand, galaxy surveys aiming to constrain PNG will require predictions that are applicable for galactic halos, whose 
masses are typically lower. Corrections to the bias due to environment and formation history may be important if non-Gaussianity is detected, 
because they could affect the value of /nl that is extracted. Our findings therefore motivate more work towards understanding the range of 
validity of J30t. and possible alte rnatives to extend the excursion set method to galactic mass-scales dD'Amico et al.l201 lahlScoccimarro et al.l 
1201 1| ; IMussq & Paraniapel2012l) . Note that we have quantified the environmental effects only due to PNG, for simplicity neglecting additional 
effects within the model that arise from non-Markovianity introduced by the filter function. Extensions of this work should focus on the full 
effects of the environment on halo formation in the context of the halo bias with non-Gaussian initial conditions. 

In addition to resolving the potential issues above, future effort should be devoted to testing the excursion set model against the bias as 
measured in N-body simulations, along with the numerous analytical predictions in the literature. We note that some authors have required 
a "fudge factor" applied to the analytical predictions in order to reduce their amplitude by ~ 10 — 25% and improve agreement with 
simulations (e.g. lGrossi et al.ll2009l : lGiannantonio & PorcianfeOld: Pillepich et al.l|2010l : IWagner & Verdell2012l) . Though this fudge factor 
has sometimes been attributed to non-spherical collapse. iMaggiore & Riotto ( 20 1 0a b i^) note th at a similar factor naturally arises within the 
excursion set model from the stochasticity of the collapse barrier (see also lRobertson e t al. 2008). We point out that the next-to-leading order 
corrections that we derived here also act to diminish the amplitude of the scale-dependent bias (i.e by < 40% for M > 1O 14 M0 /h). Though 
it is difficult to draw conclusions at this time, it is possible that the "memory" effects encapsulated in these terms, in combination with other 
effects (see below), may play a role in the reconciliation of analytical predictions with the lower amplitudes observed in simulations by some 
authors. 

In comparing the excursion set bias to N-body simulations, it is important to take into account the full range of physical ef- 
fects accessi ble to the model: 1) The effects of environment and formation history encod ed in the non-Markovianity from both non- 
Gaussianity dMaggiore & Riottoll2010"cl;[D ; Aloisio & Nataraia H ide Simone et al' I l201lj. and the current work) and the filter function 
( Maggiore & Riottcll2010al;lMa et alj|201ll). 2) Th e stochasticity of the collapse barrier to parameterize the complex nature of halo collapse 



Robertson et alj|2008t IMaggiore & Riottoll2010b]). and 3) Non-spherical collapse cha ra cterized by a depend ence of the collapse barrier on 



the filtering scale dSheth & Tormenll2002l : IZhang & Hui 2006; de Simone et al. 201 lb a; Adshead et al. 2012). With regards to an improved 
treatment of the collapse barrier, we became aware of a calculation by Adshead et al. (2012) during the preparation of this manuscript. While 
they consider only leading order terms in the expansion, (I30t . t hey add corrections to the scale-dependent bias due to a moving barrier. We 
note that they also recover the N = 3 case of equation (118) in lDesiacques et all d2011bh by setting the barrier to a constant, in agreement 
with our findings. The combination of a higher-order treatment of the connected three-point correlators with their moving barrier analysis 
would be a natural extension to this work. 
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APPENDIX A: DETAILS OF THE COLLAPSED FRACTION CALCULATION 



Since we are interested in both the leading and next-to-leading order terms, we will keep the calcu lation general for as long as p ossible, 
and then specialize to each case. We note that the approach taken here is based on the appendix of iD'Aloisio & Nataraianl d20 1 lh except 
for the following modification. Their calculation involves the collapse barrier at two different times, since they were ultimately interested in 
how primordial non-Gaussianity alters halo collapse epochs. Here, we modify the calculation to consider the collapse barrier at just a single 
epoch. The steps taken here are otherwise identical. 

We begin by separating out the terms for which at least one of i, j, and k is equal to I, 



I" i 1 1 2 2 1 3 3 

- 2~2 {S % SjS k ) c didjdk = - - 2~2 {5iSjS k )c didjdk - - ^2 {SiSjSi) c didjdi - -^2{8i8f)c didf - -{Sf) c df 



(Al) 



i,j,k — l 



l-l 



»>j=l 



- 2~2 2~2 (&<Mfc)c didjdk - i 2~2 E {5i5 3 Sk) c didjdk -^2 E (<Wj'<^}c ftftft (A2) 



i,j = l k = l+l 



i=l j,k=l+l 



i=l j = l + l 



- ^2 {SiSj6k)c didjdk - - 22 didjdi - - ^ {SiSf) c djdf . 



(A3) 



i,j,k=l+l 



i=l+l 



Consider the four terms in (IAH . For brevity, we use the following notation: Wf™ = W gln (Si; 5i+i, . . . , 8 3 ; S 3 ). Using the factorization 
property, 



W sra (8 ; 8i,...,8 s ;S s ) = W SU1 (8 ; Si,..., Si; Si) X W STa {8 v , 8 t+ i, . ..,5 S ; S s ), 



(A4) 



we find 



m-w§;: = - \wc J2 <<WA> C dtdjduwgj - i ft (wffl £ (8 1 8j8 l ) c didjWij 



i,j,k=l 

l-l 



9i (w^) E (SiSfh ftft (w*T) - \ df (w?™) £ (6i6?) c tkwfi 



\ W^d? (Wf™)) + \ d l [ft (W$?) ft (Wi 



i=l 



(A5) 



where we have obtained the first term on the right hand side by combining terms. Plugging this into equation d25t along with the contributions 
from dA2t and dA3t . and keeping only terms that are first-order in the connected 3-pt correlators yields 



P e (8 a ,S a \S l ,Si) = Ur(S c ;S l ;S a ;S s ~S l ) + P E ng (5 s ,S s \S,,Si), 



(A6) 



where we have defined 



P? e (5 s ,S s \6i,Si) 



N a + N b + N c 
nf" 1 ^; to; ft; Si -So)' 



(A7) 



with 
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(A10) 



We cannot evaluate equations JA8t through JAlOt in their current form, but instead resort to the expansion of the correlators, equation d30b 
We now consider the leading-order terms in the expansion. 



Al The leading order terms 

The leading-order contributions to N a , Nb, and N c are obtained by substituting the connected correlators evaluated at the endpoints of the 
sums in i3l\ . 
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(A13) 



We make use of identities similar to equations (48), (49), and (50) of iMaggiore & Riottol ( feOlOch (see also iMaggiore & Riottdl2010al) to 
evaluate the above expressions. Specifically, we use the following identities: 
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<9 c nf m (<5 c ; 8 ; S r , S l - S ) = ]T / d<5i . . . <Wi_i diW*™ (A14) 

-i J— oo 



1 = 1 

-1 r S 



^nf m (5 c ; S ; S t ;Si-S )=J2 dSl • ■ • d5i -! d ' d J W o7 < A15 > 

M=1 ./-oo 

d 3 c nr(Sc;S ;S l ;S,-So)= / d<Jl •• • d ^" 1 didjdhW^. (A16) 

Similarly, we would also like to use: 

s-i /•«<, 

a c nf n (<5c;<5 ! ;^;5' s -5 i )= ^ / dS l+1 . . . dd s -i cW ; s ™ (A17) 

i=!+l "'-° 

d 2 c Iir(Sc;Si;S a ;S a -Si)= / d5i+i...d5 a - 1 d i d j W!^ (A18) 

a^nf m (5 c ;^;5 a ; S. - S t ) = £ ['" Mi+i...d6.- 1 8 i d j d h Wff. (A19) 

ij,fc=i+l -°° 

Note, however, that the sum mations in equation s dA12b and dA13b are up to s and not to s — 1. The latter is required to use equations 
dATTt . <A18t . and dA19l AsrMa ggiore & Riottd d20 lOch point out, this is not a problem since we are ultimately interested in calculating 
fcoii(S s \Si, Si), which is given by 

fcoll(S s \5i,Si) = 1 - f dS s nf n (5 c ; S t ;S s ; S s - Si) - f " d8 s P e ng (J s , S s \8 h Si). (A20) 

w — oo J — oo 

We will therefore evaluate dS s iV a (j, c ) instead of N a ( bc y In that case, we can use 

d c Ue(l,s)= J2 / d5 l+1 ...dS s d t W^ (A21) 

d?U e (l,s) = / dSi+L.-dSsdidjW^ (A22) 

d 3 c U<{l,s) = ^ / d« I . H ...d*.ft0 3 -0i.W^, (A23) 

where we have defined 

US, s) = f ° d5 s nf m (<5 c ; ft; S s ; S s - ft) (A24) 

■/ — oo 

The strategy is to substitute the right-hand sides of equations dA14b - dA16b and dA21b - dA23b wherever they appear in equations dAl lb . 
dAl 2b . and dA13b . Equation dAl lb produces 

dSs at« = - ^ dl [n? m (o, i)] di [us, »)] - (sf)c did c [nf m (o, 0] a »)] - ^ d c [nf m (o, i)] 9 2 [zj e (j, *)] 

- 4^ n ? m (°- Otf s )l ~ ^r 1 # [ n " m (°' 0] 9i [us, «)] - ^ 9 i n? m (o, oaf [tf«(J, »)] . (A25) 

O 2 2 

For brevity, we have used the shorthand notation: nf m (0, 1) = Ilf m (5 c , Si — So, Si — So). Following a similar procedure for Nt and N c 
yields 

f a ds a iv 6 (1) = -Mis d i [nr(o, 0] 9c p/e(i, s)] - 9c [nr (o, 0] dl [us, s)] - (s s sf) c dtdc [n? m (o, 0] a c [us, «)] 

-<<W?)c 9c [nf m (0, 1)] did c [US, a)] 

(A26) 



© 0000 RAS, MNRAS 000, 000-000 



22 D'Aloisio, Zhang, Jeong, and Shapiro 



and 



f Sc ds s iv« = _ nf m (o, i)dl [us, a)] - nf m (o, /)d ; d c 2 »)] 

J-oo 6 2 



ftpr(o,o]a?tf«(i»»)] 
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ur\o,i)d c df[us,s)} 
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The next step is to add equations {A25}, {A26), and JA27| > and take the continuum limit. Note that if we send e — > 0, then 

U e=0 {l, s) = erf ((<5 C - 5 l )/^2(S s - &)) . 

We make use of the following properties: 



(A27) 



(A28) 



di U £ =o(l, a) = -d c U e =o(l, s) = did c U t =o(l, s) = -<9; d c U c=0 (l, s) 



(A29) 



dfU e=0 (l,s) = dlU e= o{l,s) = -did c Ut=o(l,s) 



(A30) 



diU e =o(l,s) = -d c U e=0 (l,s). 



After some rearranging, we find 
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(A31) 



(A32) 



ns,(o,o J 2 

where .A, 6, and C are defined in equations ((34) through ((36}. Equation ( IA32) may be evaluated using ( |A28t and the probability density 
llf^ (0, 1), which is given by equation i22i . Finally, performing these substitutions, differentiating, and combining with the Gaussian and 
Markovian term yields equation d33t . 



A2 The next-to-leading order terms 

Referring back to equations dA8t through JAlOt . we substitute the appropriate next-to-leading order terms from (1301 to obtain 



N, 



12) 



Ni 2) = / <Wi...<Wi_i / (Wi + i...<W,_i 



— oo 
Z-l 



d5i . . . d<5i_i / d<5 i+ i . . . dSs-i 

— oo J —oo 



gm 
I 



(1,0,0) 1-1 
1,0,0) i-1 



^(1,0,0) 

^1,1,1 



(s k - s^g^ ; 1 '] a^ (im^WT) 



I. i,i=l fc=i+l 

5 £ £ - ^)G ! ( , 1 j°; 0) + (5, + S k - 2.S. if;,;; ! ;' 1 9<(Wo7) WWfT) 

i=l 3,fc=!+l 

E £ - ^m°' 0) + (S, - S^G' , ; 1 '] ft [W^CW? 



=1 j=i+l 



d<5i . . . dft-i / " dS l+1 . . . d<5 s _i I W aj ^ (Si + Sj + S k - 3S s )d i a j d k W l a t r a 



(1,0,0) 



G 



(0,1,0) 



w 7 E (& + 5, -2&)9 i a^ 
i,j=i+i 



G 



i,i,k=i+i 

(0,0,1) 

1,1, S n2 

-o, 



wST E & - s-W, 



i=l+l 



(A33) 



(A34) 



(A35) 



The strategy is similar to our procedure in the last section. As before, we will for convenience evaluate J_ c dSsNa^ rather than the above 
expressions. However, there are some additional ingredients in the calculation due to the the presence of the Sj in the summations. We will 
encounter the quantities 
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r a = / c dSi-'-dSt-^SidiW*™, 



and 



OO t/ — oo 



r 6 = / dft/ «Mi+i...dtf,_i ftW^. 



The former may be easily evaluated by noting that 



<5c 



i = i 



r„ = ft I dS 1 ... dSi-tdiW* 1 * = ]T Si / dSi . . . ft-i W sm (ft; ft, . . . , ft-i, ft; ft) 

/ — oo 

dft +1 . . . dft_i W gm (ft; ft+i, . . . , ft; ft - Si) 



x 
j-i 



= ^ ftnf m (0; ft; ft)nf m (ft; Si; Si- ft), 

i=l 

where (see equations (79) and (80) o f lMaggiore & Riottoll2010al) 



n? m (0;ft;ft) = ^^4k ex P 

0f OS/2 



2ft 



and 



n gm (ft; ft ; ft - Si) = v^-^ 7^ JL,o ex P 



(ft - ft) 2 



2(ft - S,) 



(A36) 



(A37) 



(A38) 



(A39) 



(A40) 



Noting that IIf m (0; ft; Si)IIf m (ft; ft; Si — Si) in the sum remains finite for e — > 0, in the continuum limit we may replace the sum with an 
integral, 



dft, 



so that 



ft(ft-ft) f s < , 

r a = — l ± / ds, 

i" Jo 



s, 



ft 3/2 (ft - ft) 3 / 2 



exp 



ft 2 (ft -ft) 2 



2Si 2(S; — Si 



Using equation (108) o f lMaggiore & Riottd d201o3) . we obtain 



r /2 ft 

V 7T VS 



(2ft -ft) 2 



25, 



For T;,, we may write 



^-S s n? m (ft;ft;S 3 -S0+ ftnf m (ft;ft;S fc -SOnf m (ft;ft;ft,-S fe ) 
3 fc=i+l 



r 6 = / ds 

J — 00 

The first term evaluates to zero and, after converting the sum to an integral over dSk, the second terms gives 



ds s f Ss ds k s k - (k-W-a.) 

Js, 



(s k - ft) 3 / 2 (s s - s h fr- 



■ exp 



(ft -ft) 2 (ft -ft) 2 



2(S fc -ft) 2(S S -S*) 



which can be evaluated to 



Tfc = (ft — ft)erfc 



ft — ft 
VWs - S ; ) 



+ 



2 S; 
VS S - ft 



exp 



(ft -ft) 2 
"2(S s -ft) 



(A41) 



(A42) 



(A43) 



(A44) 



(A45) 



(A46) 
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Equations dA42t and dA46t will be needed to evaluate the final expression for the next-to-leading order contribution. 

Utilizing again the procedure of replacing sums of partial derivatives with partial derivatives of the probability density, and cumulative 
probability (i.e. equations dA14t through JA16) and ( 1A21) through dA23» . we obtain 



c d8 3 N (2) = - {diU € (i, s) [d c r a + diT a - Sid 2 c ar(o, Q - 5,9 c anr(o, Q] (A47) 

- oo 

_ dfU e (l,s) [Ta _ Sl9cUr{o l)] 
<s r G (o,o,i) 

d5 s N^ 2) = - G$° fi) d c u*(i, s) [d c r a - Sitfnr (o, i)] - -J^d 2 nr(o, [r 6 - s s d c u,(i, s)} (A48) 

-oo ' 

G (i,o,o) 

-GuT dl [( r - - StdJiriO, 0) dcU e (l, s)] - G^di [5 c n? m (0, J) (T 6 - Aft^, a))] 

* c d«5 s iv c (2) = - ^nr(o, i) (9 c 2 r, - 5* s d 3 [/ e (z, a)) - Gf; 1 ; ^, [nr(o, I) {d c r b - s s d 2 u t (i, ■))] (A49) 

G (o,o,i) 

[nr(Q) ^ _ s))] 

Plugging in r a and Ft along with Ilf™ (0, Z) and U e =o(l,s) from the last section, summing the contributes from N a , Nb and N c , and 
manipulating the result yields equation j40) - the next-to-leading order contribution to the collapsed fraction. 



APPENDIX B: DETAILS OF THE LINEAR HALO BIAS CALCULATION 
Bl Scale-independent terms 

Bl.l Leading-order 

We would like to write our final result in terms of the skewness (<Sj = {5s) c /S 2 ) rather than the three-point function, so we first convert 
derivatives of (<5g) c to derivatives of 5i 3 ' using 



dS s s dS a S, 
Let us consider the following three terms from equation ( 1591 separately 



h s ° ds s ds a {^s7~ ' (x ~ Xo) (B2) 



= &<«Hf (r^) ( " o) - (B3) 



l2 = s: {d ° )c ^s7{^sr) 



Starting with dB2t . we expand to first order in (5; about 5; = and take the limit as Si — > to obtain 



<Sl=0 



9S 3 3 V 81 S a 



Note that equation dB3b is similar to dB2t so that 



S s 3 \Si S. 

Finally, equation dB4t yields 



(B4) 



si = -s!^l ( - + i- ) 
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h = lim ( « 

s:-*oy ds t 5i=0 

Rewriting in terms of the skewness and combing 1 2 and I3, we find 



(5j)c 
2Ss 



"5? + S 8 



5? 



(B7) 



h + h = 

Combining this result with I\ gives equation l !60t . 



5 (3) 



1 . 3<5 r 2 



(B8) 



B1.2 Next-to-leading order 

Following the same procedure as in the last section, equation ( 1611 ) to first order in Si and, in the limit 5; — > 0, becomes 



G 



(1,0,0) 



f 1 

% + A 



Q^(1,0,0) „ 



We note that G^ 1 * ' can be rewritten as 



and its derivative, 



r (i,o,o) _ 1 d{S-i) c _ Si 3) f dlnSj 3) 



(B9) 



(BIO) 



dG 



(1,0,0) 



Combining these equations with dB9t gives equation ([62] 



- (3) f 91n5i 3) 

8 ^T + 4 + 



2S 2 G> 2 <sj 3) 
5i 3) 95? 



(Bll) 



B2 Scale-dependent terms 

B2. 1 Leading order 

There are four terms in equation ( 165 1 to consider: 



5S S dS a 



dS s 



(B12) 



d 2 = (S s Si) 



aj-gm /fits™ \ - 1 



dS s V c»S s 



(B13) 



6(3 - 9S S 9S S 1, 9S S J X 



(B14) 



a fgm / « fg m 

d4 = -3(tf*,) d/co11 ' /coll '° 



(B15) 



dS a V <9S S 

We begin with equation dB 13b . Expanding d-2 to first order in Si produces a constant, which will vanish when taking the cross correlation 
between 5h and 81, and a term proportional to Si. Taking the cross-correlation of Si with the term proportional to 5; gives 



(^)^(^) C (^)" 1 lirn 
where we have used the fact that 



d_ f dfZ 
dSi \ dS s 



(81 ■ 81} = (5 2 s Si} c lim (p\ 5 ,- )- lim 



dtp 



s. 



Sl=0 



(B16) 



lim i/>L -n = 



(B17) 
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to obtain the right-hand side of the equation. Conveniently, 



lim 

s.-vo 



Si 



so we are left with 



{(hSi) = {6i5i) e lim lf3\ Sl=0 



= 1, 



2S 3 



(B18) 



(B19) 



Now let us consider the contribution from (IB 121 . Again, we expand about Si = to first order in Si and take the cross-correlation with Si. 
Making use of equations ( IB 17b and JB 1 8b . we obtain 



ds 3 s^o I ds l 



dtp 



<5i=0 



Si = 



d(5 2 s 5i) c 
dS s 



(B20) 



which can be rewritten as 



{diSi) = ^r +Ss ds:{— 



(B21) 



It can be shown using the above approach that there are no contributions from ds and d&. Adding the contributions from di and d^ gives 
equation d66K 

B2.2 Next-to-leading order 
From equation ( 169) , define: 



d 5 = 2 



0,^(0,1,0) /orgm 



as 3 v 5s s J os. 



dS a 



dS s 



X 



(B22) 



(B23) 



d7 ~ aa as coil 



Si 



, _ Mo,i,o) f d f!ou,o \ d/coTi Si 
It can be shown that ds and da do not contribute. Consider the contribution from df. It is straightforward 



"Stt lim 



<5(=0 



to see that 



rgm 
Jcoll.O- 



Noting that G^f = (1/2) d(5%) c /8S 3 , 



we can re-write equation JB26) as 

as s I s s 2 dsi \ s s 



<**> = -(w) & 

Expanding dg to first order in 5; , taking the cross-correlation with Si , and the limit S s — ► 0, gives 



{daSl) = -{^sT) 



G^f lim 



which can be written as 



d 


K?) 




dSi 




s,=o_ 



s , _ (S 2 Ji)c S s d f (S 2 s Si) c 
8 ' ; 2S S 2 dS a \ S a 



(B24) 



(B25) 



(B26) 



(B27) 



(B28) 



(B29) 
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Adding the contributions from dj and ds, and substituting the three-point correlator from (|9] yields equation ( 170b . 
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